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An  Implementation  of  Adaptive  Side  Lobe  Cancellation  in 

MATLAB® 

Executive  Summary 

The  objective  of  an  adaptive  side  lobe  cancellation  system  is  to  suppress  high  duty 
cycle  and  noise-like  interference  signals  received  through  the  side  lobes  of  the  radar.  This 
is  accomplished  by  using  auxiliary  antennas  whose  antenna  pattern  approximates  to  the 
side  lobe  pattern  of  the  main  antenna.  By  suitably  phasing  the  signals  received  by  the 
auxiliaries  a  directional  anti-phase  signal  can  be  generated  that  when  added  to  the  main 
antenna  signal  “subtracts”  the  interference. 

This  report  describes  an  implementation  in  MATLAB®  of  an  adaptive  side  lobe  can¬ 
celler  system,  including  a  copy  of  the  source  code. 
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1  Introduction 

The  objective  of  an  adaptive  side  lobe  cancellation  (ASLC)  system  is  to  suppress  high 
duty  cycle  and  noise-like  interference  signals  received  through  the  side  lobes  of  the  radar. 
This  is  accomplished  by  using  auxiliary  antennas  whose  antenna  pattern  approximates  to 
the  side  lobe  pattern  of  the  main  antenna.  By  suitably  phasing  the  signals  received  by  the 
auxiliaries  a  directional  anti-phase  signal  can  be  generated  that  when  added  to  the  main 
antenna  signal  “subtracts”  the  interference. 

A  side  lobe  canceller  effectively  works  by  forming  a  beam  from  the  auxiliary  antennas 
that  points  towards  the  interference  source,  the  signal  from  this  beam  is  then  subtracted 
from  the  signal  detected  by  the  main  antenna.  The  problem  is  how  do  we  determine  the 
direction  in  which  to  point  the  beam  we  form  from  the  auxiliary  antennas. 

Mathematically  the  problem  is  to  estimate  the  complex  vector  of  N  weights  W  = 
(Wi,  W2,  ...  W /v)  to  apply  to  the  N  signals  V  =  (Vi,  V2,  ...  Vv)  from  the  auxiliary 
antennas  in  order  to  minimise  up  to  N  directional  interference  signals  in  the  side  lobes  of 
the  main  antenna  signal  (Vm).  It  may  be  shown  [1]  that  for  band  limited  noise,  appropriate 
weights  may  be  estimated  from: 


W  =  XR 

(i) 

where  M  =  E  { V*VT} 

(2) 

and  R  =  E{VmV*} 

(3) 

where  fi  is  an  arbitrary  scalar  E  { 

}  the  expected  value,  *  the  complex  conjugate, 

T  the 

transpose  and  the  inverse. 

2 

ASLC  inputs 

2.1  On  nomenclature 

The  code  nomenclature  assumes  that  the  data  has  been  Doppler  processed  and  pulse 
compressed,  but  that  is  not  essential. 

The  association  of  particular  matrix  dimensions  with  “range”  or  “Doppler”  is  one  of 
convenience  (“range”  is  the  second  dimension  and  “Doppler”  the  third  dimension  of  the 
signal  and  auxiliary  data)  and  need  not  necessarily  align  with  the  user’s  data.  However,  if 
there  is  a  change  in  the  association,  then  it  is  the  responsibility  of  the  user  to  ensure  that 
the  parameter  values  provided  are  appropriate  to  the  interpretation  of  the  data. 


2.2  Sample  data 

The  signal  data  against  which  ASLC  is  to  be  applied  is  the  Signal,  which  consists  of 
an  arbitrary  number  of  “beams”  or  “channels”  in  the  form  of  a  three  dimensional  array 
such  that  (nCs,  nRs,  nDs)=size (Signal)  where  nCs  is  the  number  of  signal  channels, 
nRs  is  the  number  of  range  bins  and  nDs  is  the  number  of  Doppler  bins. 
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The  auxiliary  data  that  is  to  provide  the  ASLC  is  Auxiliaries,  which  consists  of  an 
arbitrary  number  of  “beams”  or  “channels”  in  the  form  of  a  three  dimensional  array  such 
that  (nCa,  nRa,  nDa)=size (Auxiliaries)  where  nCa  is  the  number  of  signal  channels, 
nRa  is  the  number  of  range  bins  and  nDa  is  the  number  of  Doppler  bins.  It  is  essential 
that  nRa  ==  nRs  and  nDa  ==  nDs. 


2.3  Range  strides 

Simple  ASLC  is  implemented  with  a  single  set  of  training  data  over  the  data  set.  This 
is  not  desirable  if  the  interference  is  not  uniform  over  the  data  set.  The  baseline  code 
directly  supports  the  segmentation  of  the  data  into  “range  strides”,  with  calculation  of 
the  adaptive  weight  carried  out  for  each  stride.  The  number  of  strides  used  is  determined 
by  the  Nstride  parameter.  The  application  of  the  weights  across  the  data  is  then  by 
interpolation,  so  that  each  range  bin  has  a  unique  set  of  weights. 

In  some  circumstances  it  is  desirable  to  not  include  data  from  the  range  region  at  the 
centre  of  the  stride.  The  size  of  the  excluded  region  is  defined  by  the  RangeGap  parameter. 


2.4  Training  region 

The  extent  of  the  training  region  is  established  by  the  ClutterExtentln  vector. 

If  length(ClutterExtentIn)==4  then  the  clutter  region  is  assumed  to  extend  over 
two  regions  from  Init_bin=ClutterExtentIn(l)  to  low_bin=ClutterExtentIn(2)  and 
from  high_bin=ClutterExtentIn(3)  to  f inal_bin=ClutterExtentIn(4) . 

If  length(ClutterExtentIn)==2  then  it  is  assumed  that  Init_bin=l  and  f  inal_bin 
=nDa.  The  values  in  ClutterExtentln  are  sorted  into  an  ascending  list  so  that  they  cannot 
overlap,  although  the  regions  may  merge. 

The  sample  region  used  is  determined  by  the  SampleRegion  parameter. 

When  SampleRegion==" noise"  then  the  training  is  obtained  by  a  single  region  (from 
low_bin  to  high_bin)  in  the  data  and  the  ASLC  corrections  are  applied  over  the  entire 
data  space. 

When  SampleRegion=="clutter"  then  the  training  is  obtained  from  two  sample  re¬ 
gions  (from  Init_bin  to  low_bin  and  from  high_bin  to  final_bin)  and  the  ASLC  cor¬ 
rections  are  applied  only  over  the  training  region. 

Any  other  value  for  SampleRegion  results  in  a  selection  of  training  data  on  either  side 
of  the  middle  Doppler  bin  until  the  required  number  of  samples  has  been  obtained.  The 
resulting  ASLC  corrections  are  applied  over  the  entire  data  space. 


2.5  Training  data 

To  train  the  data,  samples  from  the  training  region  must  be  selected.  The  number  of 
samples  in  range  (Nranges  parameter)  and  Doppler  (Ndopplers  parameter)  to  be  inde¬ 
pendently  defined. 
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Radar  data  may  include  correlations  in  range  and/or  Doppler  due  to  the  waveform 
design,  oversampling,  compression  weights,  Doppler  padding  etc.  This  can  result  in  the 
use  of  correlated  data  in  the  training  which  does  not  provide  useful  training  data  (at  least 
not  in  the  baseline  implementation).  To  overcome  this  the  software  supports  the  ability 
to  require  that  only  every  nth  data  sample  is  used.  The  spacing  for  range  space  and  for 
Doppler  space  can  differ.  If  the  required  number  of  samples  will  take  the  data  outside 
of  the  region  of  validity  then  the  software  makes  appropriate  adjustments  to  the  data 
collected.  The  spacing  between  samples  is  determined  by  the  Nstep  vector.  If  Nstep  is 
a  single  value  then  the  value  is  used  for  both  range  and  Doppler,  if  it  is  an  array,  then 
the  first  value  is  used  for  the  range  dimension  and  the  second  value  used  for  the  Doppler 
dimension. 

There  is  no  checking  that  the  number  of  samples  is  sufficient  for  an  accurate  estimate 
of  the  covariance  matrix. 


2.6  Pre-computed  weights 

In  some  circumstances  it  is  desirable  to  use  pre-computed  weights  from  a  prior  data 
set.  This  is  supported  by  providing  in  the  output  from  the  function  the  weights  that  have 
been  estimated.  These  weights  may  then  be  placed  as  an  input  to  the  function  in  the 
variable  InWeights.  If  this  function  is  not  required,  then  this  variable  must  be  set  to 
InWeight s=  [] .  The  implementation  assumes  (apart  from  a  few  simple  checks)  that  the 
supplied  weights  are  appropriate  to  the  data  set  that  they  are  being  used  on. 


3  ASLC  outputs 


3.1  Output  data 

The  outputs  data  set  is  the  same  format  and  size  as  the  input  signal  data,  but  after 
ASLC  has  been  applied  to  the  signal. 


3.2  Weights  data 

Weights  data  that  is  calculated  by  the  code  is  supplied  as  an  additional  output.  This 
enables  examination  of  the  weight  magnitudes,  but  also  the  use  of  the  weights  generated 
with  a  similar  (same  size)  signal  and  auxiliary  data  (with  the  other  input  parameters  the 
same)  in  accordance  with  section  2.6. 


4  Hidden  parameters 

The  ASLC  implementation  includes  a  number  of  “hard  coded”  variables  that  may  need 
to  be  changed  for  a  particular  implementation. 
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The  parameter  ARBITARY_SCALAR  represents  //  the  arbitrary  scalar  in  equation  1  and 
is  currently  set  to  unity. 

In  training  the  data  there  is  a  simple  test  (see  source  code  lines  335-339  and  372-373) 
for  the  presence  of  a  large  target  in  the  training  data  (a  crude  cell  average  CFAR)  so  that 
large  targets  can  be  excluded  from  the  training  data.  The  magnitude  of  target  that  is 
considered  “large”  for  exclusion  is  set  by  the  value  of  LARGE_TARGET_THRESHOLD.  In  the 
code  herein  it  is  set  to  20. 

If  other  processing  has  been  applied  to  eliminate  “zero”  Doppler  effects  such  as  DPCA 
which  has  been  trained  on  the  same  data  as  will  be  used  for  ASLC,  then  it  is  (probably) 
desirable  to  exclude  that  data  from  the  ASLC  training  data.  In  the  code  herein,  the  region 
that  is  excluded  is  determined  by  the  parameter  DPCA_MASKED  and  it  ensures  that  Doppler 
bins  1 :  DPCA_MASKED  and  bins  nDa-DPCA_MASKED :  nDa  are  excluded  from  training. 


5  Limitations 

There  is  no  conditioning  applied  in  the  generation  of  the  covariance  matrix  in  this 
implementation.  Depending  upon  the  data  and  usage  of  the  code  it  is  possible  that 
the  estimate  for  the  covariance  matrix  will  be  ill-conditioned  for  inversion.  Generally 
MATLAB®  will  throw  a  warning  if  this  is  the  case. 

The  MATLAB®  code  doesn’t  make  use  of  the  (obvious)  inv(M)*R’  instead  it  uses  the 
MATLAB®  recommendation  of  M\R  ’ .  To  change  to  the  explicit  code  the  “commenting 
out”  needs  to  be  changed  in  the  source  code  at  lines  390  and  391. 

The  code  has  been  modified  to  MATLAB®  2009b.  The  majority  of  the  code  should 
work  with  prior  versions  except  where  the  dummy  variable  ~  has  been  used  in  some  return 
values.  If  the  code  is  to  be  run  on  an  older  implementation  of  MATLAB®  then  ~  will 
need  to  be  changed  to  (any)  convenient  (non-clashing)  variable  name  in  lines  125  and  129. 
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Appendix  A  MATLAB  implementation  of  ASLC 

The  following  is  the  code  used  to  compute  and  apply  the  ASLC  weights. 

1  function  [output,  weight  s]=ASLC_ARCSl  (  Signal  ,  Auxiliaries,  Nranges  ,  .  .  . 

2  Ndopplers  ,  Nstep  ,  Nstride  ,  ... 

3  SampleRegion  ,  RangeGap ,  ... 

4  InWeights,  ClutterExtentln ) 

5  %function  [output,  w  eights  ]=ASLC-ARCS1  ( Signal  ,  Auxiliaries  ,  Nranges,... 

6  %  Ndopplers  ,  Nstep  ,  Nstride  ,  ... 

7  %  SampleRegion ,  RangeGap,  ... 

8  %  InWeights,  Clutter  Extent ) 

9  % 

10  %  Compute  and  apply  the  ASLC  weights 
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%P  arameters 
%  Signal : 


35 

36 


45 

46 


Auxiliaries 
Nranges  : 
Ndopplers  : 

Nstep  : 


Nstride  : 


[Signal  Range  Doppler]  array  for  the  signal 
antenna  ( s ) 

[Aux  Range  Doppler]  array  for  the  auxiliaries 
Number  of  range  bins  to  be  used  in  covariance  estimate 
Number  of  doppler  bins  to  be  used  in  covariance 
estimate 

Number  of  range  &  Doppler  bins  between  samples  used 
in  covariance  estimate  can  be  a  vector  for  different 
step  sizes  [range  Doppler]  (not  used  if  InWeights  ~  =  [] ) 
Number  of  estimates  to  be  made  in  range  space 
SampleRegion:  Defines  the  region  within  which  training  data  is 
collected  relative  to  the  ClutterExtentln  values 
AND  the  region  to  which  the  processing  is  applied. 

If  SampleRegion=="  noise  11  then  the  training  data  is 
collected  in  the  region  low-bin :  high-bin  (see 
ClutterExtentln)  and  is  applied  to  the  ENTIRE  Doppler 
extent  of  the  data.  If  SampleRegion=="  Clutter  ''  then 
training  data  is  collected  in  two  regions 
Init-bin  :  Low-bin  and  high_bin  :  final-bin  (see 
ClutterExtentln)  and  the  processing  is  only  applied  to 
the  same  region.  Other  values  for  SampleRegion  define 
a  default  space  in  the  middle  of  the  Doppler  coverage, 
extent  of  the  range  gap  in  around  the  stride  central 
range  that  is  to  be  used  in  estimating  the  weights 
(not  used  if  InWeights  ~  =  [] ) 

Precomputed  weights  ,  if  weights  are  to  be  estimated 
from  the  data  then  this  MUST  be  a  null  (=[])  array, 
if  it  is  NOT  null  ,  then  the  value  will  be  used 
(subject  to  some  sanity  checks)  overwrite  Nstride 
and  as  the  ASLC  weights.  The  form  of  InWeights  is 
assumed  to  match  the  form  of  the  output  parameter 
weights  . 

ClutterExtentln  A  Either  a  four  value  array  [ Init-bin  low-bin 

high_bin  final-bin]  that  defines  the  clutter  region 
to  be  [  Init-bin  :  low-bin  ]  and  [high_bin  :  final-bin  ]  , 
or  a  two  value  array  [low-bin  high_bin ]  that 
assumes  Init -bin=l  and  final-bin^nD  where  nD  is  the 


RangeGap : 


In  Weights  : 
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50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 


%  maximum  Doppler  bin.  Note  the  final  list  is  always 

%  sorted  into  assending  order. 

%  output:  [Signal  Range  Doppler]  array  after  ASLC  processing 

%  weights:  [stride  Aux]  array  of  the  weights. 

% 

%  Notes  : 

%  1.  Code  syntax  assumes  post— Doppler  processing . 

%  2.  There  is  no  conditioning  of  the  covariance  matrix  prior  to 
%  inversion  in  this  code!  MATLAB  may  throw  warning  messages ! 

% 

%  Original  by  APShaw,  2f  July  2009 

%  Modified  for  clutter  processing  (RANGE-GAP)  by  APShaw,  10  Aug  2009 

%  Modified  for  pre— computed  weights  in  clutter  by  APShaw,  21  Aug  2009 

%  Modified  to  exclude  very  edge  Doppler's  from  '' edge  11  option 
%  by  APShaw,  01  Sep  2009 

%  Modified  to  speed  up  pro  cessing  by  eliminating  squeeze  functions  in 
%  the  weight  computations  using  a  mixture  of  a  squeeze  and  a  reshape. 
%  by  APShaw,  08  Sep  2009 

%  Modified  to  add  ''  clutter  11  specific  option  by  APShaw,  08  Sep  2009 

%  Clean  up  and  improve  comments  to  form  version  for  ARCS  study 
%  distribution  and  set  the  coding  conform  to  MATLAB  2009b.  Eliminate 
%  un—used  SampleRegion  options.  Replace  fixed  offsets  with 
%  parameterised  values  (DPCAJMASKED  and  LARGE-TARGET-THRESHOLD) . 

%  Changed  '  Clutter  Extent  11  variable  to  enable  the  clutter  extent 
%  region  to  be  tricked  into  other  regions  or  to  limited  regions  of 
%  the  clutter.  Reformated  to  suit  incorporation  into  tech  note. 

%  by  APShaw,  15  Jun  2010 

% 

% 

% 

%  (c)  Copyright  ,  Commonwealth  of  Australia  ,  2009,  2010 


%% - 

%  Fixed  parameters  and  early  initialisations 
% - 


%  arbitrary  scalar  applied  to  noise  cancellation 
ARBITARYJdCALAR-  1 ; 

%  Set  the  threshold  to  be  used  in  determining  if  there  is  a  target 
%  present  in  the  training  data  that  needs  to  be  excised.  Value  is 
%  the  SNR  (linear) 

LARGK  TARGKT  TIIRESI  IOI.I)= 2  0 ; 

%  Set  the  number  of  bins  to  be  eliminated  from  the  Doppler  region  due 
%  to  the  effects  of  DPCA 
DRCA  MASKI-’I)  3; 

%  Gap  in  range  data  to  be  used  at  center  of  each  stride  ,  to  prevent 
%  training  on  ’’self  data”. 

RANGE_GAP=RangeGap ; 
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103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 


%  Identify  the  number  of  signal  and  the  range/Doppler  map  size 
[nS  nRs  nDs]=  size  (  Signal  ) ; 

%  Identify  the  number  of  auxiliaries 
[  nA  nRa  nDa]  =  size(  Auxiliaries  ); 

%  Initialise  the  output  array 

output=zeros  (size  (Signal))  +  lj*zeros  (size  (Signal)); 


%% - 

%  Validate  the  input  data 
% - 


%  determine  if  we  are  using  pre— computed  weights  and  validate  the  data 
%  provided . 

if  ~  isempty  ( InWeights  ) 

%  check  that  the  number  of  weight  channels  equals  the  number  of  aux 
%  channels 

if  ndims  ( InWeights)==3 

[nWs  ,  nWc]=  size  ( InWeights  ) ; 
assert  (and  (nWs  =  nS  ,  nWc  =  nA ) 

'Number^of -weight  ^channels  „must~match„SIG~and„AUX„channels  '  ) ; 

else 

[~  ,  nWc]=size  (InWeights  ) ; 
assert  (nWc=nA ,  .  .  . 

'Number^of  ^weight  ^channels  JVlUSThmatch ^number u of  WIBCchannels  ' ) 
nWs  =1; 
if  nS>l 

%  if  there  are  more  signal  channels  than  weights  assume 
%  applies  only  to  the  first  channel. 

Signal  ( 2  :  end  , :  ,:)  =  []; 

end 

end 

Nstride=nWs; 

COMPUIEWEIGinS=false  ; 

else 

COMPUIEWHGmS=t  rue; 

end 


%  Check  number  of  signal/ auxiliary  Doppler  s  and  ranges  match 
assert  (  and  (nRa  =  nRs,  nDa  =  nDs)  ,  ... 

'  Range  „and  ^  D  opp  ler  ^  bin  „size„of„signal  „and  ^auxiliaries  JVIUSThmatch  '  ) ; 


%  Check  the  Clutter  Extent  values 
if  ~ isempty  (  ClutterExtentln  ) 

if  length  (  C  lu  t  ter  Extent  In  )==2 

ClutterExtent  =[1  ClutterExtentln  ( 1 )  ClutterExtendln  (2)  nDs]; 
elseif  length  (  ClutterExtent In)>=4 
C  lut  ter  Ext  ent=Clut  ter  Extent  In  ; 

else 
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156  disp(  '  Clutter  „ extent  „b  ad  ly  u  cl  efined^— ^using^default^region  '  ); 

157  SampleRegion= '  default  '  ; 

158  ClutterExtent  =[1  nDs  nDs  nDs  ]  ; 

159  end 

160  ClutterExtent=sort  (  ClutterExtent  ) ; 

161  else 

162  disp(  'Clutter„extent~is~not~defined~—  ~u sin g„default~ region  '  ); 

163  SampleRegion= '  default  '  ; 

164  ClutterExtent  =[1  nDs  nDs  nDs]; 

165  end 

166 

167  %  Set  the  values  for  the  step  sizes  in  range  and  Doppler  depending 

168  %  upon  the  value  (s)  passed  by  Nstep. 

169  if  length  ( Nstep)==l 

170  DNranges=max(  Nstep  ,  1 ) ; 

171  DNdopplers=max(  Nstep  ,  1 ) ; 

172  else 

173  DNranges=max(  Nstep  ( 1 )  ,  1 ) ; 

174  DNdopplers=max(  Nstep  ( 2 )  ,1); 

175  end 

176  DNranges_lialf=floor  (  Nranges  .  *  DNranges  .  /  2  ) ; 

177  DNdopplers_half=  floor  (Ndopplers.*  DNdopplers  .  /  2 ) ; 

178 

179 

180  %  Check  the  number  of  ranges/ dopplers  is  within  the  size  of  the  arrays 

181  assert  (  (( Nranges  .*  DNranges))<=  nRa  ,  .  .  . 

182  '  V alue^of  u<Nranges*  Nstep>„  larger  ..than  „number„  range  „bins  ,  i nv.;dat  a,  '  ) ; 

183  assert  (  Ndopplers  .*  DNdopplers<=  nDa  . 

184  '  V  alue  „of  „<  Ndopplers  *Nstep>„  larger  „than„doppler„bins„in„data  '  ) ; 

185 

186  %% - 

187  %  Determine  regions  to  be  used  for  estimations 

188  % - 

189 

190  %  compute  the  center  range  bin  index  for  each  stride 

191  start=DNranges_half  +  l; 

192  stop=max(nRa— DNranges_half  ,  start  ) ; 

193  check  =  [] ; 

194  if  Nstride>l 

195  strides=floor(start+(  stop  — start  )*(0:Nstride—  l)/(Nstride  —  1 ) ) ; 

196  %  remove  any  strides  with  the  same  value  (to  prevent  redundant 

197  %  calculations  and  interpolation  failing) 

198  check_count  =0; 

199  for  count=2:  Nstride 

200  if  strides  (count)==strides  (count— 1) 

201  clreck_count=check_count  +1; 

202  check  (  check_count)  =  count  ;  %#ok<AGROW> 

203  end 

204  end 

205  if  check_count 

206  strides  ( check  )  =  [] ; 

207 

208  else 
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strides=floor  (nRa .  /  2  ) ; 

end 

Nstride=length  (  strides  ) ; 

%  compute  the  centres  of  the  Doppler  samples  and  generate  a  list  of  the 
%  Doppler  bins  to  be  used  for  weight  estimation . 

Doppler_step=max(  1  ,  DNdopplers  ) ; 

SampleRegion=lower  (  SampleRegion  ) ; 
switch  SampleRegion 
case  'clutter  ' 

%  Make  a  list  equal  to  the  entire  Doppler  space 
Doppler  _listO=l:  Doppler  _step  :  nDs  ; 

%  Eliminate  those  bins  that  are  NOT  in  the  clutter  region 
Doppler _  1  i s 1 0  (... 

and  (  Doppler_listO  >  ClutterExtent  (2)  ,  ... 

Doppler_listO  <  ClutterExtent  ( 3 )))  =  [] ; 

Doppler_listO  (  Doppler_listO  <  ClutterExtent  ( 1  ))  =  [] ; 
Doppler.listO  (  Doppler_listO  >  ClutterExtent  ( 4 ))  =  [] ; 

%  Use  the  eliminated  region  as  the  Doppler  list  for  training 
Doppler_list=Doppler_listO  ; 

%  remove  the  DPCA_masked  region  from  the  training  list 
Doppler_list  (  D o p p  1  e r _  1  i s t  <=DPCAJVL4SKED )  =  [] ; 

Doppler.list  (Doppler  _  1  i  s  t  >=nDs— DPCA  MASKED )  =  [  ] ; 
case  'noise  ' 

%  Make  a  list  equal  to  the  entire  Doppler  space 
Doppler  _listO=l:  Doppler  _step  :  nDs  ; 

%  Eliminate  those  bins  that  are  in  the  clutter  region 
Doppler  .list  0  (  Doppler  .list  0  <=  ClutterExtent  ( 2 ) )  =  [] ; 

Doppler  _  1  i s 1 0  (  Doppler.listO  >=  ClutterExtent  (3))  =  [] ; 

%  Use  the  eliminated  region  as  the  Doppler  list  for  training 
Doppler_list=Doppler_listO  ; 
otherwise  %'  default  ' 

D0=floor  (nDs  .  /  2  ) ; 

Doppler.list  =  ... 

max(  1  ,D0—  DNdopplers.half  ) :  ... 

Doppler_step  :  min(DO+DNdopplers_half  ,nDa) ; 

end 

NumberDopplers=length  (  Doppler  _  1  i s t  ) ; 

%  Failsafe  provision  —  ensure  we  have  at  least  one  Doppler  in  the  list 
if  NumberDopplers<l 

Doppler_list=floor  (nDa .  /  2  ) ; 

NumberDopplers  =  l; 

end 


%% - 

%  Estimate  the  weight  sets  for  each  stride 
% - - - - - 

if  COMPUIEWEIGmB 

%  initialise  the  weights  and  the  intermediate  matrices 
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weights=zeros  (  Nstride  ,nA,nS)  +  lj  .*  zeros  (  Nstride  ,nA,nS); 

M_init=zeros  (nA)  +  l  j  .  *  zeros  (nA ) ; 

R_init=zeros  (nS  ,nA)  +  lj  .*  zeros  (nS  ,nA) ; 
tot al_power=ones  (nS  ,  1 ) ; 

%  loop  over  the  number  of  range  strides  to  take  over  the  data  set 
for  stride_counter  =1:  Nstride 

%  Determine  the  extent  of  region  where  range  data  is  to  be 
%  taken  ,  making  accomodation  for  the  gap  in  the  range  data 
%  and  retaining  the  specified  number  of  range  bins  in  the 
%  data  set . 

Rangel=st rides  (  st r ide_count er )—  DNranges_half— RANGE_GAP; 

Range2=  strides  (stride_counter)  +  DNra n ge s _h al f+RANGE_GAP ; 

%  If  the  range  extents  go  beyond  the  limits  of  valid  data, 

%  adjust  them  to  be  within  the  range  limits 
if  Range2>nRa 

Excess=Range2— nRa ; 

Rangel=Rangel— Excess  ; 

Range2=nRa ; 

end 

if  Rangel <1 

Excess=l— Rangel ; 

Range2=Range2+Excess  ; 

Rangel  =  1; 

end 

%  Double  check  the  effects  of  the  offsets  incase  a  large 
%  number  of  range  bins  are  in  use. 

Range  l=max  ( 1  ,  Range  1 ) ; 

Range2=min  ( Range2  ,  nRa ) ; 

%  Determine  the  extent  of  the  region  where  there  is  a  gap 
%  in  the  range  data. 

Skipl=strides  (  stride_counter )— RANGE_GAP; 

Skip2=st rides  (  stride_counter )+RANGE_GAP; 

%  Check  the  guard  region  doesn't  go  beyond  the  limits  of  valid 
%  data  and  if  they  do  adjust  them  accordingly 
if  Skip2>nRa 

Excess=Skip2— nRa ; 

Skipl=Skipl— Excess  ; 

Skip2=nRa; 

end 

if  Skipl<l 

Excess=l— Skipl  ; 

Skip2=Skip2+Excess  ; 

Skipl  =1; 

end 

%  Double  check  in  case  a  large  value  of  RANGE-GAP  is  in  use. 
Skipl=max(  1  ,  Skipl  ) ; 

Skip2=min  ( Skip2  ,  nRa); 
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%  Make  a  list  of  the  range  bins  to  be  used  for  this  stride  and 
%  then  remove  those  that  are  in  the  skip  region 
ValidRanges=Rangel  :  DNranges  :  Range2  ; 

ValidRanges  ( and  ( ValidRanges>=Skipl  ,  ValidRanges<=Skip2  ))  =  []; 

%  Fail-safe  provision:  ensure  that  we  have  at  least  one  range 
%  bin  left  in  the  stride. 
if  length  ( ValidRanges)<l 

ValidRanges=st rides  (stride_counter  ); 

end 


%  generate  a  vector  of  total  signal  powers  in  this  stride 
%  accomodating  the  gap. 
for  counter  =  l:nS 

total _power  (counter  ,:)  =  squeeze  (  .  .  . 

sum  ( sum  (  Signal  (nS  ,  ValidRanges  ,Doppler_list)  ... 

.  *  conj  (  Signal  (nS  ,  ValidRanges  ,Doppler_list  ))))); 

end 

cells  =squeeze((  length  (ValidRanges  ).*  length  (  Doppler.list ))  —1); 

%  initialise  weight  computation 
M=M_init  ; 

R=R_init  ; 

%  compute  the  covariance  matrix  for  this  range  stride 
for  range  _counter=l:  length  (  ValidRanges  ) 
r ange_bin=V alid Ranges  (range.counter  ) ; 

%  Speed  up  functionality  :  extract  data  for  this  range  bin 
T_Signal=squeeze  (Signal  (:  ,  range  _b  in  ,  :  )  )  ; 
if  nS==l 

T_Signal=T_Signal  .  !; 

end 

T  _  Auxiliar  ies=squeeze  (Auxiliaries  (:  ,  range  _b  in  ,  :  )  )  ; 

%  speed  up  functionality  :  reshape  to  avoid  ’’squeeze”  inside 
%  the  inner  loop  . 

t  .sign al=resh ape  (  T -Signal  ,  size  (  T -Signal  ,  1 )  ... 

*  size  (  T-Signal  ,2)  ,1); 
t  _  Auxiliar  ies=resh  ape  (T.Auxiliaries  ,... 

size  (  T-Auxiliaries  ,  1 )  *  size  (  T.Auxiliaries  ,2)  ,1); 

for  doppler_counter  =1 :  Number  Dopplers 

doppler_bin=D  op  pier -list  (doppler_counter  ) ; 

Svector=t  .signal  ((doppler.bin—  l)*nS  +  l:  .  .  . 

(  doppler.bin  —  l)*nS+nS  ) ; 

Spower=Svector  .  *  conj  (  Svector  ) ; 

%  If  there  isn  't  a  large  target  present  in  the  test 
%  cell  then  add  this  cell 

i  f  ~  (max(  (  Spower  .  /  ( ( tot  al  .power —Spower  )./  cells  )... 
>LARGE_TARGET_THRESHOLD ) ) ) 


11 


DSTO-TN-0955 


368  Avector=t -Auxiliaries  ((  doppler_bin —l)*nA+l:  ... 

369  (  doppler_bin  —  l)*nA+nA) ; 

370  M=M  +  (  Avector )  *  Avector  '  ; 

371  R=R  +  (  S  vector  )*  conj  ( Avector  ).' ; 

372  %  count er=counter  + 1  ; 

373  end 

374  end 

375  end 

376 

377  %N ormalis ations  turned  off  to  save  computational  cycles 

378  %M=M.  /  counter ; 

379  %R=R . /  counter  ; 

380 

381  %  Compute  the  weights  :  note  MATLAB  recommendation  to  not 

382  %  use  inv  function. 

383  %  weights  (  stride-counter  inv  (M)  *R' ; 

384  weights  (  stride_counter  ,:)=M\R'; 

385  end 

386  %  need  to  take  complex  conjug ate  of  the  weights  —  this  is  because 

387  %  we  have  allowed  MATLAB  to  do  the  conjug  ate  transpose  instead  of 

388  %  the  transpose  as  this  is  marginally  faster  on  the  test  machine. 

389  weights=conj  (  weights  ) ; 

390  else 

391  %  If  we  are  not  computing  the  weights  from  the  data  ,  then  use  the 

392  %  weights  that  have  been  supplied  in  the  input  parameters . 

393  weights=InWeights ; 

394  end 

395 

396 

397  %% - 

398  %  Apply  the  weighted  auxiliaries  to  the  signal 

399  % - 

400 

401  %  Range  values,  using  ’’start”  and  ’’stop”  (calculated  previously )  values 

402  %  to  avoid  extrapolation 

403  ranges=min  (max  ( 1 :  nRs  ,  start  )  ,  stop  ) ; 

404  %  the  Doppler  expansion  of  the  range  weights  to  simplify  the  weight 

405  %  application 

406  expander=ones  (nDs  ,  1 )  ; 

407  %  initialise  the  correction  matrix 

408  correctionO=zeros  (nRs  ,  nDs)  +  lj  *  zeros  (nRs  ,  nDs  ) ; 

409  %  initialise  the  output  matrix 

410  %output=zeros  ( nS ,  nRs ,  nDs)  +  lj  *  zeros  (nS ,  nRs  ,  nDs ) ; 

411 

412  %  Apply  the  weights  to  each  signal  in  turn 

413  for  counterO  =l:nS 

414 

415  %extract  the  weights  for  the  current  signal 

416  Sweiglrts=squeeze  (  weights  ( :  , :  ,  counterO  ) ) ; 

417 

418  %  initialise  the  correction  for  this  signal  channel 

419  correction=correctionO  ; 

420 
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%compute  the  weight  to  be  applied  at  each  range  bin  ( output  weights 
%array  has  dimensions  nA  x  nRs 
if  Nstride==l 

%  one  weight  applies  to  every  range  bin 
Wout  =  ((  squeeze  (Sweights(l  ,:))).  ' )  *  ones  ( 1  ,  nRs  ) ; 

%  compute  the  correction  factor 
for  counter  =  l:nA 

correction=correction+  ... 

squeeze(Auxiliaries(counter  ,:)).*  ... 

( (  expander  *Wout(  counter  ,  : ) )  .  '  )  ; 

end 

else 

%  Different  weights  for  each  range  bin 
for  counter  =l:nA 

%  interpolate  the  computed  strides  to  each  range  bin  for 
%  each  array  for  this  signal  channel 
Win=squeeze  (  Sweights  ( :  ,  counter  ) ) ; 

Wout  (  counter  ,:)=  inter  p  1  (strides  ,  Win,  ranges  ,  '  spline  '  ) ; 

%  compute  the  correction  factor 
correction=correction+  ... 

squeeze(Auxiliaries(counter  ,:)).*  ... 

( (  expander  *Wout(  counter  ,  : )  )  .  '  )  ; 

end 

end 

%  apply  the  correction  to  the  sigal ,  including  the  arbitrary  scalar 
%  adjustment  of  the  weighted  correction . 
switch  SampleRegion 
case  ' clutter  ' 

%  Initially  make  the  output  equal  to  the  input 
output ( counterO  , :  , : )  =  squeeze ( Signal ( counterO  , :  , : ) )  ; 

%  Apply  the  correction  only  over  the  entire  clutter  region 
%  Note  —  this  does  include  the  DPCAJAASKED  region. 
output  (  counterO  , :  ,  Doppler _listO  )  =  . . . 

squeeze(Signal(counterO  ,  Doppler.listO))—  ... 
correction  (:  ,  Doppler_listO)  *  ARBIT ARY JSC  ALAR ; 
otherwise 

output ( counterO  , :  , : )  =  squeeze ( Signal ( counterO  , :  , : ) )  —  . . . 
correction  .  * ARBITARY_SCALAR ; 

end 

end 


13 


DSTO-TN-0955 


14 


Page  classification:  UNCLASSIFIED 


DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION  1.  caveat/privacy  marking 
DOCUMENT  CONTROL  DATA 

2.  TITLE  3.  SECURITY  CLASSIFICATION 

An  Implementation  of  Adaptive  Side  Lobe  Can-  Document  (U) 

cellation  in  MATLAB®  Title  (U) 

Abstract  (U) 

4.  AUTHORS  5.  CORPORATE  AUTHOR 

A.  P.  Shaw  Defence  Science  and  Technology  Organisation 

PO  Box  1500 

Edinburgh,  South  Australia  5111,  Australia 

6a.  DSTO  NUMBER.  6b.  AR  NUMBER  6c.  TYPE  OF  REPORT  I  7.  DOCUMENT  DATE 


DSTO-TN-0955  014-078  Technical  Note  July,  2010 


8.  FILE  NUMBER 

9.  TASK  NUMBER 

10.  SPONSOR 

11.  No  OF  PAGES 

12.  No  OF  REFS 

2010/1076861/1 

07/044 

PM  C  &  W 

14 

1 

13.  URL  OF  ELECTRONIC  VERSION  14.  RELEASE  AUTHORITY 

http://www.dsto.defence.gov.au/corporate/  Chief,  Electronic  Warfare  and  Radar  Division 

reports/DSTO-TN-0955.pdf 


15.  SECONDARY  RELEASE  STATEMENT  OF  THIS  DOCUMENT 

Approved  For  Public  Release 

OVERSEAS  ENQUIRIES  OUTSIDE  STATED  LIMITATIONS  SHOULD  BE  REFERRED  THROUGH  DOCUMENT  EXCHANGE,  PO  BOX  1500, 
EDINBURGH,  SOUTH  AUSTRALIA  5111 _ 

16.  DELIBERATE  ANNOUNCEMENT 

No  Limitations 

17.  CITATION  IN  OTHER  DOCUMENTS 

No  Limitations 

18.  DSTO  RESEARCH  LIBRARY  THESAURUS 

Radar,  Radar  interference,  Radar  jamming, 

Signal  Processing,  Adaptive  signal  processing, 

Side  lobes, 

Software,  Computer  programs 

19.  ABSTRACT 

This  report  describes  an  implementation  in  MATLAB®  of  an  adaptive  side  lobe  canceller  system, 
including  a  copy  of  the  source  code. 

Page  classification:  UNCLASSIFIED 


